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We discuss the three neutrino flavor evolution problem with general, flavor-diagonal, matter po- 
tentials and a fully parameterized mixing matrix that includes CP violation, and derive expressions 
for the eigenvalues, mixing angles and phases. We demonstrate that, in the limit that the mu and 
tau potentials are equal, the eigenvalues and matter mixing angles Q\2 and #13 are independent of 
the CP phase, although #23 does have CP dependence. Since we are interested in developing a 
framework that can be used for S matrix calculations of neutrino flavor transformation, it is use- 
ful to work in a basis that contains only off-diagonal entries in the Hamiltonian. We derive the 
"non-adiabaticity" parameters that appear in the Hamiltonian in this basis. We then introduce the 
neutrino S matrix, derive its evolution equation and the integral solution. We find that this new 
Hamiltonian, and therefore the S matrix, in the limit that the mu and tau neutrino potentials are 
the same, is independent of both 623 and the CP violating phase. In this limit, any CP violation 
in the flavor basis can only be introduced via the rotation matrices, and so effects which derive 
from the CP phase are then straightforward to determine. We show explicitly that the electron 
neutrino and electron antineutrino survival probability is independent of the CP phase in this limit. 
Conversely, if the CP phase is nonzero and mu and tau matter potentials are not equal, then the 
electron neutrino survival probability cannot be independent of the CP phase. 

PACS numbers: 14.60.Pq 



'Electronic address: Kneller0ipno.in2p3.fr 
^Electronic address: Gail_McLaughlinancsu.edu 



2 



I. INTRODUCTION 

Neutrino flavor transformation in vacuum [1-4] and in matter [5, 6] continues to be the subject of much attention 
both experimentally and theoretically. On the experimental side we have moved from a situation where little was 
known about the mixing parameters roughly a decade ago to one where half of the mixing parameters are reasonably 
well known. The experimental status upon the mixing parameters is that 5m{ 2 = &.0t_o.i eV 2 , 1.9 x 1CT 3 cV 2 < 
\dm% 3 \ < 3.0 x 1CT 3 eV 2 , sin 2 20 12 = 0.861^, sin 2 2<9 23 > 0.92, and sin 2 26 13 < 0.19 (at 90%) [7, 8]. The sign of 
<5m 2 3 , whether # 2 3 is less than or greater than 45°, and the extent of CP violation or the Majorana phases arc not 
currently known. 

On the theoretical side the focus while #i 2 was completely unknown was primarily upon solar neutrinos but after 
the results of SNO [9] development of neutrino flavor transformation theory has increasingly focused upon supernova 
neutrinos which display much richer phenomena. Two primary shifts in supernova neutrino flavor transformation 
theory have occurred. First, simplified and static profiles have increasingly been abandoned, in favor of more realistic, 
dynamic and turbulent density profiles [10-19]. Secondly, the implications of more complete descriptions of the 
potentials are now being studied, such as neutrino self interactions [20-51] and the "matter mu-tau" potential [42, 44, 
53]. These studies have uncovered new behaviors of the neutrinos, and in the process in has become apparent that 
the calculational tools employed in the past have limited applicability to these newer situations. 

While some progress has been made in joining together these two new developments [54], in this work we are 
interested primarily in the first. An approximation that is often made in treating the evolution of the neutrino 
wavefunction through more realistic density profiles is that the mixing of the neutrinos can be treated in a scries 
of two flavor mixing schemes even though there are three neutrino flavors. The computationally desirable reduction 
from three to two flavors is based upon the observation that typically only two out of the three neutrino states will be 
mixing with each other at a resonance and that resonances are widely separated. But with the use of more realistic 
density profiles multiple, closely spaced, resonances can occur and it has been shown [14, 16, 19] that finely grained 
energy dependent phase effects appear in two flavor mixing. This suggests that the phenomenon also extends to three 
flavors, although this has not yet been investigated. These effects can only be seen when the phase information of 
the wavefunctions is retained throughout the entire calculation, so most approximation methods will miss them. The 
full range of effects of the CP phase [55-57] can also be hard to determine with approximation methods and these 
effects arc purely three-flavor phenomenon. Finally, the small difference between the matter potentials for the [i and 
r flavors may also become important even in cases where the neutrino potentials are small, because this potential 
difference produces a third matter resonance. 

As an alternative to the usual integration of the neutrino wave functions, we consider the S matrix formulation of 
neutrino flavor transformation which was shown to be computationally more efficient in the two flavor case [16]. In 
this paper, we derive a convenient basis which can be used with the 5* matrix, which we call the "adiabatic" basis. 
We shall consider the full, generalized three flavor mixing including all phases and all three matter potentials. In 
section §11, we derive the eigenvalues and discuss their behavior. In section §111 we present expressions for the matter 
mixing angles and phases and the Hamiltonian in the matter basis. We then derive the three flavor Hamiltonian in 
the "adiabatic basis" [16] - the basis where the Hamiltonian is completely off-diagonal - in section §IV and from it we 
discover the expressions for the three-flavor non-adiabaticity parameters. We demonstrate how the S matrix is found 
from this Hamiltonian in section §V, its equation of motion and integral solution and discuss its parameterization. We 
then discuss the two flavor approximation and observe some identities of neutrino propagation that have implications 
for the CP phase. Finally in section §VI we present our conclusions. 



II. THE FLAVOR BASIS 

In this section we write down the Hamiltonian which describes neutrino flavor transformation in the flavor basis 
without approximations and determine its eigenvalues. 

The 3x3 Hamiltonian in the flavor basis is composed of two terms: the rotated vacuum Hamiltonian U KU^ 
and the matter induced potentials V^'s which are diagonal in the flavor basis if we ignore neutrino collective effects 
[58]. The full Hamiltonian is 



HU) = UKU^ + V {f) 



(1) 
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where 

/ ki \ j /m? \ 

K = =E+— mi , (2) 

\ fc 3 / 2A \ m 2 3 J 

/V e \ 

=0^0, (3) 
\0 V T J 

and i? is the neutrino energy, mi the neutrino masses and V a the matter potentials that are, possibly, functions of 
position. We have included here three matter potentials for the neutrinos for generality. We could have removed one 
of these potentials, since we are free to add to the Hamiltonian an arbitrary multiple of the unit matrix including a 
term that is a function of position since the only effect of such a term is to introduce a phase. As a consequence of 
this property only the relative difference between potentials (and eigenvalues) are important for observable quantities, 
not the absolute values. However, we choose not to impose this feature from the start, and instead allow it to 
appear automatically as we proceed. For neutrinos, the difference in the potentials V e and is the well known 
SV efl w V2 Gf (n e - — n e + ) [6] where Gf is Fermi's constant and n e is the electron number density. For antineutrinos, 
the potentials change sign. The potential difference between and V T arises due to radiative corrections to neutral 
current scattering and is smaller than SV efJ , and 5V eT by a factor of <~ 10~ 5 in typically encountered matter. Thus, 
the splitting is usually only important for neutrino propagation through supernova profiles at high density. 

In equation (1) U is a unitary matrix that relates the flavor and mass bases. The mixing matrix U has nine elements 
but the unitary conditions mean that four elements may be expressed in terms of the remaining five after specifying 
the phase of the determinant. The unitary conditions also place two further restrictions upon the magnitudes of the 
five independent elements by establishing two relationships between them. Thus, in general, U is parameterized by 
three magnitudes and six phases. (See e.g. Ref. [59] for discussion of the construction of an SU(3) matrix.) It is 
traditional to select the independent elements to be those in the first row and last column but other choices are also 
valid. The three independent magnitudes, which must all be smaller than unity, may be expressed in terms of three 
mixing angles 6*i 2 , #13 and 6*23- The form for the matrix U that we use for neutrinos is 



/l \ / C12C13 si2Ci3 aia \ (e~ mi \ 

U = e l5 -si 2 c 23 e l£ - ci2Si 3 s 2 3 c 12 c 23 e te - s 12 s 13 s 23 c 13 s 23 e - *" 2 (4) 

\ e< f3+5 '> ) \ s 12 s 23 e le - c 12 s 13 c 23 -c 12 S2 3 e li - s 12 s 13 c 23 c 13 c 23 J \0 e""* 3 / 

where C12 = cos6>i 2 , S12 = sin#i 2 etc. and ct\, a 2 , a 3 , [3, 5 and e are the six phases. For antineutrinos, instead of U 
we take U* for the unitary transformation matrix. Other than the three a's we have tried to make the expression for 
U as similar as possible to that found in [7] but note that in this more general parametrization the CP phase is not 
S, but is instead e. The phases /3 and S may be absorbed by redefinitions of the neutrino fields in the standard model 
Lagrangian showing that their absolute values are not observable. This degree of freedom can be used to recover the 
usual convention for CP by setting 5——e. Even though the values of these two phases cannot affect observables we 
will continue to keep them because, for consistency, during the calculation we must keep track of their derivatives. 
The mixing matrix U in equation (4) can be written as the product of six terms 



U = B(p,S)e 23 (e2 3 )E(e)e 13 (6 13 )Q 1 2(e 1 2)A(a 1 ,a2,a 3 ), 
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e lS 
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(5a) 



(5b) 
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With this unitary matrix we find that the elements of H^ arc, in full, 



H ee = V e + (kic\ 2 + k 2 s\ 2 ) c 2 3 + k 3 s\ 3 , 

H eli = e~ lS [(h ~ fc i c i2 - k 2s\ 2 ) C13S13S23 - (fci - k 2 ) c 12 s 12 c 13 c 23 e~ le ] , 



[(k 3 - kic 2 12 - k 2 s\ 2 ) C13S13C23 + (fci - k 2 ) c 12 s 12 c 13 s 23 e l£ ] , 



= Vp + (fci + k 2 ) C23 + k 3 cf 3 sl 3 + (fcic 2 2 + k 2 s\ 2 ) (s\ 3 sl 3 - cj 3 ) + 2 (fci - k 2 ) c 12 s 12 s 13 c 23 s 23 c e , 
H^ T = e~ tf3 [{k 3 - k x s\ 2 - k 2 c\ 2 ) c 23 s 23 - (k 3 - k X (? Y2 - k 2 s\ 2 ) s\ 3 c 23 s 23 + (fci - k 2 ) c 12 si 2 s 13 (c 23 e te 



s 23 e 



V T + (fci + fc 2 ) s 23 + k 3 c 13 c 23 + (fcic 12 + k 2 s 12 ) (s 13 c 23 - s 23 ) - 2 (fci - fc 2 ) ci 2 s 12 Si 3 c 23 s 23 c e . 



(6a) 
(6b) 
(6c) 
(6d) 

% 
(6e) 

(6f) 



The remaining elements are given by H^ e = H*^, H Te = H* T and H Tfl — H* T but, since is diagonal, all the 
off-diagonal elements are equal to the vaccum values i.e. H e ^ = H e ^ etc. Note also that ct\, a 2 and a 3 do not appear 
in these expressions, although the phases <5, (3 and the CP phase e do appear. We also point out that + H TT , 
H^ — |-ff^r| 2 and \H ell \ 2 + \H eT \ 2 are all independent of e. 
The characteristic polynomial for the eigenvalues fci , fc 2 and fc3 of H can be written as 



kf-fk 2 + I + 3Q j ~h- ( % + + 2Rj =0 



(7) 



where ki is any one of the three eigenvalues and T, Q and R are three functions equal to 



Q = — 

^ 18 



(^H ee — Hfifj^J + (^H ee — H TT ^j + (jl^fi — H TT 



\H e J 2 + \H eT \ 2 + \H, < 2 



(8a) 
(8b) 



R — ^ (2 H TT — H ee - (2 - H ee - H TT ^j (2 H ee — H^ - H TT ^j + - {H efl H^ T H re + H eT H Tfl H^ e ) 



\H^ T \ 2 (2H ee — H uu — H TT ) + \H ( 1 
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) + \H eT \ 2 (2 H^ — H ee - H TT ^j + \H efJ _\ 2 (2 H TT — H ee — H Mt ^j 



(8c) 



All three functions T(V e , V^, V T ), Q{V ei V^, V T ) and R(V e , V^, V T ) are polynomials in V e , V^, and V T and T, Q and R 
at a given V e , and V T are independent of the basis. In full these polynomials are 

f = T + Ve + Vp + Vr, (9a) 

Q = Q + 2 q e V e + 2 q„ V„ + 2 q T V T - ^ [dV^ + 5V 2 T + SV 2 e ] , (9b) 
R = R + r e V e + rpVp + r T V T -q e (V 2 + 2V f ,Vr)-qp(V 2 + 2V e V T )-q T (V 2 + 2V e Vp) 

+^(^ + SV eT )(6V^ e + SV^ T )(6V Te + SV TV ), (9c) 

where T, Q and R are the values of T, Q and R in the vacuum, SV a p = V a — Vp, SV 2 p = (V a — Vg) 2 and the six 
functions q e , q^, q T and r e , r M , r T are 



Qe = 

^ = 

(h = 

r e = 

r u = 



1 
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{H^n + H TT — 2 H ee ) , 



(H ee + H TT — 2 H^p) , 
(Hee + Hpfj, — 2 H TT ) , 

(H 

(Hpp - H TT ) 2 + (H^n - H ee ) 2 — 2 (H TT - H ee ) 2 + 3 (\H ^ T \ 2 + \H^ e \ 2 — 2 |-ffre| 2 ) 

t- \(H TT - Hee) 2 + (H TT — H^^) 2 — 2 (Hee — H^^) 2 + 3 ( \H Te \ 2 + \H TI1 \ 2 - 2 |H eA j| 2 ^j . 



18 

1 

18 

1 

18 
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H^) 2 + (Hee - H TT ) 2 - 2 (H^ - H TT f + 3 ( \H e ^\ 2 + \H eT \ 2 - 2 |fl- MT | 2 



(10a) 
(10b) 
(10c) 
(lOd) 
(lOe) 
(lOf) 
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In the above expressions, we use H ee etc. to refer to the first element of = U K W, i.e the flavor basis Hamiltonian 
in vacuum. Since the quantities, T, Q and R are independent of the basis and their vacuum values are most easily 
expressed in the diagonal, mass basis. In this basis their values are 




where Skij = ki — kj and Skjj = (ki — kj) 2 . The discriminant D = Q 3 + R 2 is negative definite consequently all three 
eigenvalues are real. In terms of T, Q and R the eigenvalues k\ , ki and k% arc [60] : 



appropriate value in vacuum, i.e ki becomes ki in the vacuum and similarly k 2 — > k 2 and k% — > ^3. As the orderings 
of the vacuum values depend on the hierarchy, there are two possible choices of values for these angles u\ , lu 2 and 
CJ3 . In the case of the normal hierarchy (NH) where k\ < k 2 < k%, the angles must be u>\ = 2n, u 2 = 4-7T and 0J3 = 0, 
while in the case of the inverted hierarchy (IH) where k\ < k\ < k 2 , the angles must be u\ = A-k, uo 2 = and W3 = 2n. 
Antineutrinos have the same vacuum eigenvalues as neutrinos, so these angles are the same for both neutrinos and 
antineutrinos. 

In the limit where V^=V T =0 we find that the characteristic polynomial, equation (7), is independent of the CP 
phase. It is clear that the vaccum values of T, Q, and R, lla-llc, are independent of the CP phase, but a careful 
inspection of equations (9a) - (9c) using previously established identites between the elements of the Hamiltonian 
shows that the e dependence drops out here as well. This implies the eigenvalues are also independent of e. 

In addition to the hierarchy, the behavior of the matter eigenstates is strongly influenced by whether #23 > 7r/4 
or #23 < tt/4. Thus, when considering all the matter potentials, our present knowledge of the neutrino mixing 
parameters allows four possibilities for the general evolution of k\, ki and £3 [42] with V. These are shown in figure 
(1). The solid vertical line in the middle of each panel indicates vacuum with V e — 0, so the left portion of each panel 
corresponds to antineutrinos and the right portion to neutrinos. Note that in the figure, the axes are scaled linearly, 
not logarithmically. To make the figure we set both 6*12 and #13 to be small and included a non-zero potential for V T 
proportional to, but smaller in magnitude than, V e and set = 0. In reality, #12 ~ 34 degrees, but the behavior of 
the eigenstates is still qualitatively similar to what is shown in the figure. Also is not zero, but in order to best 
illustrate the general behavior of the eigenvalues in a figure, we have used the freedom of adding to the Hamiltonian 
a multiple of the unit matrix. 

When the mixing angles are small each eigenvalue tends to track one of the potentials. At "resonances" they switch 
which potential they follow therefore ensuring that the ordering of the eigenvalues is invariant. With three potentials 
there are three resonances commonly referred to as the "L" , "H" and the "/it" . The L resonance, L standing for "low 
density" , occurs at the smallest positive value of V e and, because Sm 21 is known to be positive, the two states that mix 
are always the neutrino states V\ and v 2 . The H resonance, H standing for "high density" occurs at the intermediate 
value of |V e |. Which two states mix depends upon the hierarchy: for a normal hierarchy (NH) it is states i> 2 and £3, 
for an inverted hierarchy (IH) it is the antineutrinos states v\ and P3. The two states that mix at the [it resonance, 
which occurs at the highest value of |V e | depends both upon the hierarchy and whether #23 < tt/4 or #23 > tt/4. For 
a NH and #23 < tt/4 it is states v\ and v 2 but if #23 > tt/4 then it is states i>i and ^3, in an IH states v\ and i>i mix 
when 023 < tt/4 but v\ and £3 when #23 > t/4. Out of the four possibilities the only one which causes two eigenvalues 
to mix twice, when considering only a monotonically decreasing density profile, is the scheme of a NH and #23 < t/4 
because one sees that states i>\ and v-i mix once at the fj,r resonance and again at the L resonance. This double mixing 
of v\ and v 2 raises the possibility that interference effects might appear even in monotonic profiles. We will examine 
the consequences of this in future work, but now we use the results of this section to compute the Hamiltonian in the 
matter basis for use with the S matrix. 




(12) 




The three angles u>i, 0J2 and 0*3 are chosen to ensure that each eigenvalue takes on its 



III. THE MATTER BASIS 



In this section, while keeping track of all the angles, phases and their derivatives, we write down the Hamiltonian 
in the matter basis. We then make a few remarks about the large density limit. 
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FIG. 1: The four possibilities for ordering and evolution of the eigenvalues fei, &2 and £3. The potential V M is set to zero while 
V T oc V e but \V e \ > \V T \. The solid vertical line is located at V e = 0. The vacuum mixing angles 812 and #13 were set to small 
values so that the eigenvalues closely followed the potentials. In all four panels the solid line is fci, the short dashed line is ki 
and the long dashed line is £3. In the top two panels the hierarchy is normal, and in the bottom pair the hierarchy is inverted 
(IH). For the left pair of panels the vacuum mixing angle #23 obeys #23 < 7r/4, for the right pair #23 > 7r/4. The units of either 
axis are linear but arbitrary. 



With the eigenvalues determined we can begin by deducing the angles #12, #13, #23 and phases (3, S and e in the 
unitary transformation U that relates to the matrix where k±, ki and ^3 appear on the diagonal i.e 

(k x \ 

fj\ H (f) (J = K= i ~k 2 . (13) 

\o hj 

After making the change of basis we find that the Schrodinger equation is 
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^ -(*-<)*<-'. 

= ff( m ) V (m) . (14b) 

The term W dll / dx appears because the eigenvalues are functions of position which requires that U also be a function 
of position. Obviously if W dU /dx = then the matter basis Hamiltonian would be diagonal but this can only occur 
if the density is constant otherwise, in general, ij( m ) has non-zero off-diagonal elements that come from dU/dx. 
In order to evaluate H^ m > we need to compute W dU /dx. 

The matter mixing matrix U is parameterized in exactly the same way as U in equation (4) i.e. as 



1 °^ \ /ci 2 ci3 Si 2 5i 3 _ s 13 e \ / e-^ 

U = ( e lS -S12C23 e" - C12S13S23 C12C23 e le - S12S13S23 C13S23 e~ la2 |(15) 

e l (P+V J V S i2*23e J£ - Ci 2 s 13 C23 -Ci 2 s 2 3 e 4e - S12S13C23 ci 3 c 23 / \0 e^" 

where Cy = cos% etc. After some rather lengthy algebra we find that the matter mixing angles 6*i 2 , 6\ 3 and 6*23, 
and the matter phases /3, 6 and e can be expressed in terms of the eigenvalues and the elements of the flavor basis 
Hamiltonian as 



tan 2 6*12 



<5fci 3 




-k)(Hrr 


-fc 2 )-|H M r| 2 


Sk 2 3 




- k)(Hrr 


"fcl)-|^r| 2 " 



(16a) 



• 22 (H^-h)(H TT -h)-\H^ T \ 2 

sin z 6»i 3 = — — -~ — -~ - — , (16b) 

ok 13 dk 23 

tan 2 9 — ( ^ Te — ^H±i^Il ~ J [ ^ er — H efJ _(H TT — k 3 ) \ (16c) 
\-ff T p — H Te (H[j_n — k 3 ) J \H efl H^r - H eT (H^^ — k 3 ) J 

- H Te H e u — Hu e H eT — (H TIJ/ — H„ T )(H ee — k 3 ) ,.„,> 

1 tan/5 = - - - = ^- L , (16d) 

H Te H efl + H^ e H eT — (H TI1 + H^ T ) (H ee — k 3 ) 

1 Hn T H Te — H eT H T n — (Hu e — H e n)(H TT — k 3 ) 
ztano = — - — = ^- L , (16e) 

H^t H Te + H eT H T ^ — (H^ e + H e)Jl ) (H TT — k 3 ) 
. _ <5fcl2^fcl3<5fc23Cl2Sl2C 2 3 Si3C23S23 . ,, rc \ 

sine = — = = sine. (161) 

5fci2l5/Cl3<5fc23Cl25l2C 1 3Sl3C23S23 

These equations match those already written down in the literature in the appropriate limits. The last equation, (16f), 
for e is the Naumov [61] and Harrison & Scott [62] identity and comes from the fact that the off-diagonal elements 
of H given in equations (6b), (6c) and (6e) are not functions of position. The Toshev [63] and Kimura, Takamura & 
Yokomakura [64] identities are no longer valid when SV^r ^ 0. Other than a re-arrangement of the ordering of the 
eigenvalues and some notational changes the expressions for tan 2 #12 and sin 2 9\ 3 are the same as those one may find 
in Bellandi et al. [65] ; our expression for tan 2 #23 differs from Bellandi et al. because we allow the possibility of e 7^ 
and do not set (3 = 0. When performing numerical calculations, the expressions (16a) - (16f) become increasingly 
difficult to use when the densities become large. In section III A we discuss this limit, how one can derive asymptotic 
expressions for equations (16a) to (16f) and demonstrate why the numerical difficulties arise. 

The most straightforward way to use these equations is to first derive the eigenvalues ki,k 2 , k 3 at each density and 
then determine the matter phases and angles. But since these are trigonometric expressions, more than one value of 
the angles or phases can satisfy the identity. The matter angles 612 and 6*13, always remain the quadrant in which they 
are first defined, which removes any ambiguity in using equations (16a) and (16b). Furthermore, these two angles 
are independent of eiiV^ = V T . This can be observed because we have already established that the eigenvalues and 
the expressions H^ + H TT , H^ H TT — \H^ T \ 2 are independent of e. However 6*23 can migrate with V into the second 
quadrant if its vacuum value is defined to lie in the first quadrant so we do have an ambiguity to resolve when using 
equation (16c) to find 23 . An examiniation of equation (16c) shows that 6*23 = 7r/2 when 

£3 = %^ = fc 9 o, h = H„ ^JL^L = fcgV (17) 
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In general fcgo is complex and so at no density can we obtain an eigenvalue £3 that satisfies the condition and, therefore, 
9 2 3 never wanders out from the first quadrant. But when H e/1 H Te is pure real (which requires the CP violating 
phase e to be zero or 7r) then ^90 is also pure real and therefore there is a density at which #23 approaches and then 
passes 7r/2. Thus we see that the evolution of #23 with the potential depends upon e in a way that #12 and 613 do not. 
If we have a case where e is either zero or tt then by keeping track of the behavior of these denominators, equation 
(17), we can resolve the ambiguity and determine the correct quadrant for 823- If we have set #23 in the vacuum to be 
less than tt/2 and fc 3 < k 90 then #23 must become greater than tt/2 for k 3 > k g0 - If on the other hand we have set #23 
in the vacuum to be less than tt/2 but fc 3 > fcgo, then §23 must remain in the first quadrant for all positive definite V. 
Similar arguments can be made for the zeros of the numerator of equation (16c) that become relevant in the inverted 
hierarchy. One finds that the numerators vanish at the root fco and, again, fc is generally complex except when e is 
either zero or tt. So, if #23 in the vacuum is greater than zero and fc 3 > k then #23 < if fc 3 < fco- The CP phase 
e can range from zero to tt, but since we can determine not only sine but also cose there is no ambiguity in the CP 
phase. In principle, there are ambiguities in the phases (3 and 5, but in practice, only their derivatives are used not 
their absolute values. Finally we note there are no expressions for Si, 612 and <S 3 so we are free to pick anything for 
them including functions of x. Finally, once the matter mixing angles and phases are found we note that, like U, the 
mixing matrix U can be written as the product of six terms 



U 



Bifi, 5) 623(^23) E(e) 013(013) 612(612) A(a ly a 2 , a 3 ), 
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(18a) 



(18b) 



Next we need to calculate the derivatives of these angles and phases. A surprisingly simple result is that the 
derivatives of the eigenvalues are given by the expression 



dki 
dx 



2 

dx 



2 dV^ 
dx 



+ 



2 dV T 

dx 



(19) 



By unitarity then dk\/dx + dk2/dx + dks/dx = 
the invariancc of the trace of the Hamiltonian. 



dV e /dx + dV^/dx + dV T /dx which is also an identity that stems from 
With this result we find 
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(5fci3<5fc23 

Sk 12 



C23S23 Ski 3 Sk 2 3 

Cl2Sl2C23Se 5&12 



2 dSV^ 2 

S23 ^x~ + C2 " 

_ 2 d8V elJ , 



23 dx 
d6V f 



T <-23 



dx 



dx 



C23S23 
Ski 3 Sk23 



Cl2Sl2C 1 3C23S23CeSk 12 



dSV a 



dx 



(c 2 2 5k 13 + s 2 2 6k 23 ^J 



dx 



'23" 



~2 



dx 



^ 2 d5V eT 

+ ° 23 ^x~ 



S13S23 <5fci3<5fc 2 3 
Cl2Sl 2 Sl3S£ 5k 12 



n d5V eu . o d8V e ., 



dx 



C23S23 Ski3Sk 2 3 



~ 2 d5V efl 
323 dx 



7.2 



^ 23 ^7 
d6V„ 



5 2 - 

A 23 



dSVur 



dx 



+ - 



C23S23S £ 



dx 

12 



C12S12S13 i5fci2<5fci 3 (5fc23 



{^12^ 12 Ci 3 Skf 2 ~ Ci 2 Sl 3 <5fcl3 Sl 2 sf 3 Sk 23 ^ 



dSV. 



11T 



dx 



(20a) 
, (20b) 
(20c) 
(20d) 
(20e) 

(20f) 



Note that d/3/dx, dS/dx and de/dx are all proportional to sine and it is the appearance of this term in equation (20f) 
that validates our previous statement that e = for all V if e = 0. One also sees that when 5V^ T = there is a 
very simple relationship between df3/dx, dS/dx, de/dx and d6 23 /dx the last of which can be easily solved to give the 
Toshev identity [63]. 
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So we have finally reached the point where we are able to show that the matrix W dU / dx has the form 

„~ T ,2 /0 e lS& ^ 0\ ,2 /0 c 12 e lS& ^ 

m dU = dfa[ _ e _ lS&i2 Q o I +^ I S 12 e^ 23 



dx \ J dx \ -C12 e^ 5 " 13 -§i2 e" 4 " 5 " 23 

-2ici 2 si 2 si3Se Sis (5 2 2 e 4? + «i 2 e~' g ) e 4<5 " 12 -S12C13 e^" 13 ' 

+ I - S ^(c 2 12 e-^ + rs 2 12 e tf -)e-^ 2id 12 § 12 § 13 S e c 12 c 13 e 45523 -» 



S12C13 e -* 55l3+l? -C12C13 e- t<5523+l? 



d/3 .+ - d5 



1 - cf 2 cf 3 -Cl2Sl 2 Cf 3 e 



12^13 






-ci 2 ci 3 si 3 e " l5ai3 

-5i 2 Si 2 e^ 12 
-ci 2 s 12 e 15512 c 2 2 

oy u:t y / y / v o o 

(21) 

where the vector up is given by 

up = ( (si 2 s 23 e" - 5i 2 si 3 c 23 ) e~ iai , - (ci 2 s 23 e" + si 2 si 3 c 23 ) e~"* 2 , ci 3 c 23 e" l<i3 ) . (22) 

So far <Si, o;2 and «3 are unconstrained which permits us to choose anything for them. It is at this point that we make 
our selection that di, a 2 and 6:3 be constants. Since when using the scattering matrix approach, it is convenient to 
work with a Hamilitonian that has zeros along the diagonal, another, sensible, alternative would be to assign di, d2 
and d 3 so that the diagonal elements of dU /dx were exactly zero i.e. the contributions from da±/dx, da 2 /dx and 
dotz/dx in the equation above exactly canceled the contributions from d0 2 s/dx, d/3/dx, dS/dx and de/dx that appear 
along the diagonal. But, as we shall shortly show in sec IV, our desire to cancel the diagonal elements of W dU /dx 
can be achieved more conveniently by another route hence we opt to set the three phases to be constants. 



A. The High Density Limit 

In the limit of high density, the expressions for the matter angles and phases, (16a) to (16f), become difficult to 
work with. The on-diagonal components of the flavor space Hamiltonian, H ee , and H TT , are dominated by the 
matter potentials which are much larger than the off-diagonal components. Further, the eigenvalues k\, k 2l and k% are 
also dominated by the matter potentials so when using the expressions (16a) to (16f) one is sometimes in the situation 
where two large numbers must be subtracted to find a small difference. For this reason it is useful to find expansions 
for the eigenvalues in powers of the potentials. In this section we derive the asymptotic limits of the neutrino and 
antineutrino eigenvalues when V e — > 00 assuming V e 3> and V e ^> V T . 

We begin by considering the two quantities Q and R. The original expressions for these quantities may be rewritten 

as 



Q = 



—V 

v e 

9 

-V 2 



R = 



9 

V 3 

■ e 

27 



(lSge + Vy + Vr) K - W + v? - 9 [Q + 2 g„y„ + 2 g T V T }) 

Ve + V 2 

1 + ^ + -^ 

Ve V 2 

_ 3(18 ge + ^ + K) _ 3 (V 2 - 4V^V T + V 2 - 18 [r e - 2 g^V T - 2 g r V^}) 
2 V P 2 K 2 



(23) 
(24) 



+ 



{2V 3 - W 2 V T - 3V^V 2 + 2V 3 + 54 [R + +r T V T - 2q e V„V T - q^V 2 - q T V 2 ] ) ' 



2 V 3 



V 3 
27 



1 + 



3qq 
2V R 



V} 



V 3 



(25) 
(26) 
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For antineutrinos 



Q 



-v 2 

9 

-V 2 



1 , (18 q e -V,- Vr) { (V? VyVr + V? 9 [Q 2 q,V, 2 g T y T ]) 



9 

R = — ls- 
27 



v; T v 2 



U 2 



(27) 
(28) 



1 | 3(18g e -y M -V;) 3 (V 2 - AV^Vr + V 2 - 18 [r e + 2 g M K + 2 



214 



2 V 2 



(2V 3 - 3V M 2 V T - 3V^V 2 + 2V T 3 -U[R- r„V M - r T V T - 2q e V^V T - q^V 2 - q T V 2 ] ) 



2 V? 



V 3 

_ e 

27 



3a Q 6r Cfl 
2V e + V} V} 



(29) 
(30) 



Using these expressions we find that cos a; — ► 1, i.e. u — > in the V e — > oo limit. If we expand out coso; in terms of 
powers of u) then 



UJ 2 LU 4 

1 -T + 24 + 



1 + a R /V e + b R /V 2 + c R /V e 3 



(l + a Q /V e + b Q /V 2 ) 



2VV2 ' 



(31) 



Solving for u and expanding out the right-hand side of the resulting expression up to order 1/V e 2 we are eventually 
led to 



oj = 



3a 2 Q + 12b Q -8b R [ 5 a 3 Q + 12 a Q {b Q -b R ) + 8 c R ] 



2 V, 



AV 2 ^3a 2 Q + 12b Q -8b R 



Ve V 2 



which defines the quantities and b u . When we substitute in the definitions of clq, 6q etc. we find that 
3 V3 



a u = 



bu> = 



4 H, 



\ 2 +(H^-H TT y 



(32a) 
(32b) 

(33a) 



3\/3 r- 



+ Htt — 2H e 



A\H^ T \ 2 +(H^-H rT y 



3V3 


2 {H e ^H ^ T H re 


+ HerHr^H^e) + ( \H e ^ 2 — \H eT \ 2 




2 




\Ja\h»t\ 2 + (h^-h tt ) 2 





(33b) 



For antineutrinos we find instead that cosw — > — 1 which indicates the angle u approaches n as V e — > oo. In this 
case we instead expand cos w in terms of cD — n and find the asymptotic limit is 



v vr 



(34) 



where a w and b u are given by exactly the same expressions as and b u , although the elements of the Hamiltonian 
are, in this case, taken from the Hamiltonian appropriate for antineutrinos. With the asymptotic expression for u) 
and u> determined we can then proceed to the eigenvalues. For neutrinos these are given by 



r (T + Ve + Vy + Vr) , 2V e 
1 = 3 "3" C ° S 



OJ + UJi 



^/l + a Q /V e + b Q /V 2 . 



(35) 



Using the expression for the expansion of oj then, up to 1/14, we eventually find 



11 



A V e ( LUi — TT 

— — COS 



T + Vp + Vr a Q /wa 2a u . (u t \ 

H £ h -rr cos — — sin — 

3 3 V 3 / 9 V 3 / 

3 6q + y ] cos (^) + (a Q a w + 2 6 W ) sin (^) 



1 

"9V e 



AV e 



COS 



4 

-~>i - 7T 



COS 



3 

CJi + TT 



Ci 



V, 



(36a) 
(36b) 



Notice that if u>i — 2tt or Ui = An then the term linear in V e vanishes: it only survives when u>i = 0. Substituting in 
for a w etc. we find 



Cbi = 



1 + 2 cos 



u, 



1 — COS 



(I) 



A\H^ T \ 2 + (H^-H TT f 



V3 



sin (|), (37a) 



2 (H elx H^ T H Te + HerHr^H^e) + (\H efJ _\ 2 — \H eT \ 2 ) — H TT ^j 



+ 



(\H e ,\ 2 + \H eT \ 2 ) cos(^) 



(37b) 



Note that when Ui = then Ci = iJ ee so that the first two terms of this eigenvalue h are h = V e + H ee = H e 
The expression for the eigenvalues in the case of antineutrinos is instead 



= 3—^ + — cos 



U + LOi 



l + a Q /V e + b Q /V 2 . 



(38) 



When we substitute in the asymptotic expansion for u we find 



r 4 V e f (di - TT 

ki = cos 

3 



6 



cos 



Ui + 3ir\ (T-V^-Vr) a Q fui + n\ 2a 



AV e 



1 

9Ve 
COS 



cos 



sin 



Ui + TT 



3 ~ a Q or a*\ 

4 ~ Q + Y) cos 



(aqa w + 2 bj) sin 



Ui + TT 



Ui - TT 



COS 



Ui + 3 TT 

6 



V, 



(39a) 
(39b) 



Here the linear term in V e vanishes if Ui 7^ 2tt. In this case the expression Ci is equal to 



Q = 



1 — 2 cos 



Ui + TT 



1 + cos 



Ui + TT 



{a 



>A\H^\ 2 + (H^-H TT ) 2 

+ - = sin 

V3 



Ui + TT 



«t = 



2 (H e ^H^ T H Te + H eT H Tfl H^ e ) + (\H efJ _\ 2 - \H eT \ 2 ) - H TT ^j 



sin 



(\H e ^\ 2 + \H eT \ 2 ) cos 



Ui+TT S 



(40a) 



(40b) 



We can now write down the eigenvalues in the high density limit. For neutrinos and the normal hierarchy u\ — 2tt, 



12 



LU2 = 47r and W3 = so in the infinite density limit 



«3 



+ 



^r| 2 + 



+ ^ + o(i/K 2 ), 



i/ ee + ^ + o(i/K 2 ). 



For antineutrinos in the normal hierarchy 

61 



^3 



fcl 



ff ee + ^ + o(i/v e % 



1 



+ 



N 





2 + 




— i/ TT ^ 




4 




2 + 




— H TT ^j 




4 


r, ^2 = 


= and w 3 


= 2tt 



+ ^ + o(i/K 2 ), 



+ ^ + o(i/K 2 ). 



+ 



\H UT \ 2 + 



(^Hfifi — H T rj 



+ ^ + o(i/K 2 ), 



k 2 - #ee + TT + 0(1/K 2 ), 



+ g + o(i/K 2 ). 



And for antineutrinos in the inverted hierarchy 



^3 



03 



»,.x| 2 + 



+ ^ + o(i/K 2 ), 



+ 



\ 



»,.x| 2 + 



X - 

L + ^ + o^), 



//., +^ + o(i/y e 2 ). 



(41a) 

(41b) 
(41c) 

(42a) 
(42b) 
(42c) 

(43a) 
(43b) 

(43c) 

(44a) 

(44b) 
(44c) 



The expressions we have derived for k\, ki and fc 3 in the asymptotic limit can be employed together with equations 
(16a) to (16f) to more easily determine the correct values for the mixing angles and phases in situations where the 
matter potential is large. 



IV. THE ADIABATIC BASIS 



Returning to the problem of neutrino propagation we now introduce a new basis %l>^ a \ which we call the adiabatic 
basis, with = W{x)4> {a \ 

After making this change of basis we find the Schrodingcr equation has become 



dx \ dx dx 

(45a) 

iJ< a V (a) - (45b) 
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and we choose W so that it removes the diagonal of H^ a \ The matrix W is simply 

/ exp[-2z7n/>i] \ 

W= exp[-2i7r0 2 ] (46) 

\0 exp[-2z7T0 3 ] J 

where 4>\ , </> 2 and 4>3 are defined to be 

#i 
dx 

#2 

dx 

#3 
dx 

If we had used di, d 2 and d 3 to cancel off the diagonal elements of W dU /dx then equations (47a) to (47c) would 
only contain the first terms in the above expression, i.e. only the contribution proportional to fci, fc 2 and k 3 . The 
name "adiabatic" comes from the consideration of the two-flavor oscillation problem where the phase analogous to 
the 0i's is the adiabatic phase of the matter states. For three flavors since their are multiple phases which also contain 
additional terms proportional to dSV^/dx the situation is more complex. However, this basis is useful for providing 
insight into the phase effects, e.g. Ref [16], and the effects of the CP phase e as described in Sec. VB. Further, it is 
an ideal basis in which to perform neutrino flavor transformation calculations using the S matrix prescription. After 
removing the diagonal elements of we can write out as 




Sl 2 Sl 3 C 23 S 23 S e Sk 23 dSV. 



5i2 <5fci 2 5A:i3 dx J ' 

Cl2Sl3C 2 3§23Se $kl3 ddV^r \ 

§12 6ki2 Sk 2 3 dx J ' 

Cl2gl25i 3 C 2 3g 23 S e Ski2 dSV^r 

s~i3 Sk 13 Sk 23 dx 



(47a) 
(47b) 
(47c) 



ij(o.) = 



Skl2 p g 2 J7T(5012 l Sk 13 p g 2 V7T64>13 



6 k 







2.7T * 13 1 

5/C23 p g2nrS<j>23 



(48) 



\ Z 7T 



13 1 



13 „ Sk 2 3 p* e -2in54>23 
2 7T 23 



where, as usual, Sfaj = 4>i — (j>j- This equation defines three functions Ti 2 , Ti 3 and T 23 which are the non-adiabaticity 
parameters for 3 flavor neutrino oscillations. By matching the expressions we find 
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2?re 



Sk x 



ddl2 +si3(^ + si 2 e-^ d ° 23 



dx 



dx 



23 



dS 



de 



-i (si 2 s 23 e le - £i 2 si 3 c 23 ) (ci 2 s 23 e 4£ + si 2 si 3 c 23 ) — - zc^s^c^ ^ - «Ci 2 Si 2 — 



ri 3 = - 



2 Tresis 




<W 13 


<5fci 3 


\Cl2 


dx 


2 7re lS&23 




d0 13 




I Sl2 


Sk 2 3 


dx 



- si 2 ci 3 e 



d9 23 . t - - - % i - - - s - - dp _ _ _ 
+ I (si 2 s 23 e - c i2 si 3 c 23 J c i3 c 23 , — - ic 12 c 13 s 13 



dx 



dS_ 

dx 



- - -a d9 23 \ _ _ dp _ _ _ dS 

+ c\2C 13 e — i (ci 2 s 23 e + si 2 si 3 c 23 J c 13 c 23 , — - «si 2 ci 3 si 3 — 



dx 



dx 



(49a) 
(49b) 
(49c) 



Note how the derivatives of the two phases j3 and S appear in these expressions. When we substitute in the expressions 
for the derivatives we find 



2ne 



12 



13 



Sk\ 2 



- • i 2 dSV etl | . 2 dSV e . 



ci 2 si 2 c 13 s 23 



dx 



(c 12 s 12 (c 23 - s 23 ) + si 3 c 23 s 23 (c 12 e le - s 12 e te )) 



dSV, 



dx 



2ne 



5k\ 3 



Cl 2 Ci 3 Si 3 S 23 



~2 dSV efl . 2 dSV e , 

+ c 23" 



dx 



dx 



Sl2Cl3C 2 3S23e 



r23 — 



2ire 



tSa 2 3 



5k 23 



Sl 2 Cl 3 Sl 3 s 23 



_ dSV e ^ . 2 d5V e , 

+ c. 



dx 



23" 



dx 



ci2Ci 3 c 23 s 23 e 



dx 
,;d5V UT 



dx 



(50a) 
(50b) 

(50c) 
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If we restict ourselves to 8V^ T = then we see that Ti 2 oc C\ 2 si 2 / 8k\ 2 dV e / dx which is exactly the two flavor 
non-adiabaticity parameter as described in [16]. The other similar terms, Ti 3 oc ci 3 Si 3 /Skf 3 dV e /dx and T 2 3 oc 
Ci3Si3/<5fc| 3 dV e /dx also take on similar meanings. This Hamiltonian is constructed in such a way that, similar to the 
two flavor case, the focus is on the "non-adiabatic" pieces of the solution, i.e. the places near the resonances where 
the matter eigenstates are likely to swap. There are some corrections for the case of three flavors that come from the 
terms cf 3 , c\ 2 and §12 and the three, arbitrary, complex exponentials in ri2,Ti3 and 1^23 respectively. It can be seen 
that it is the behavior of 612 that selects between mixing of states v\ and D 3 or between D 2 and £3. 

V. THE SCATTERING MATRIX 

The Schrodinger equation for the evolution of the neutrino wavefunction is 

^ = -iH(x)i>(x). (51) 



When we integrate equation (51) we obtain 



4>{X) = V(*o) - * / d Xl H{ Xl ) V(zi). (52) 
We choose the initial point to be at Xq. Repeated substitution of this result into itself yields 

%I>{X) = i/>(X )-% ( dx 1 H 1 ^(X ) + (-i) 2 f d Xl H x { 1 dx 2 H 2 4>{X ) + . . . , (53a) 

J Xq J Xq J Xq 



= 1 1 - 1 f dx l H l + {-if ( d Xl Hi f 1 dx 2 H 2 + . . . \ %I>{X Q ) 

I J Xq J Xq J Xq J 



(53b) 



where the subscripts on the H's mean Hi = H(xi). This equation defines the scattering matrix S(X,X ) since 

i>{X) = S(X,X )i ! (X ). (54) 

The upper limits on the integrals appearing in equation (53b) indicate the space ordering but we can change all the 
upper limits to X by using identities such as 

f d Xl Hi f 1 dx 2 H 2 = ^- f thu f dx 2 {#i H 2 Q{ Xl - x 2 ) + H 2 H 1 Q(x 2 - n)} (55) 

•J Xq J Xq ' J Xq J Xq 

where ®{x\ — X2) is the Heaviside step function. Using this result and similar identities for the higher order multiple 
integrals, allows us to write S as 

r X 1 \2 r x r x 



S(X,X ) = l + (-i) f dxiH! + LA. [ dxi f dx 2 T(H 1 H 2 ) 

J Xq ' J Xq J Xq 

+ f dx i I dx 2 f dx 3 T(H 1 H 2 H 3 ) + ..., (56) 

v Xq J Xq J Xq 

where T is the space/time-ordering operator. Now that it is defined we can simply insert equation (54) into the 
Schrodinger equation and find that S also obeys the differential equation 

i^=HS (57) 
ax 

which describes 9 coupled equations for the elements of S. From this equation we can also derive that the phase of 
the determinant \S\ — e J * is simply. 

$(X, X ) = - [ Tt(H) dx. (58) 
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So we find that in the adiabatic basis & a '(X, X ) = and does not vary because the Hamiltonian is traceless. Since 
the initial condition is that S^(X a ,X ) = 1 we see that in the adiabatic basis S is a member of SU(3). For all bases 
S has the property that it obeys the product rule 

S(X, X Q ) = S(X, X,)S(X„ X ). (59) 
and the probability that a neutrino with initial state \Dj(X )} is later detected as state \i>i{X)) is 

P(\i> j )^\i>i)) = \S ij (X,Xo)\ 2 . (60) 
S is, in general, a member of U(3) and this restriction means that any element of S(X,X ) satisfies the relationship 

Sij(X, X ) = e^ x ' x <>> X ) (61) 

where C 4 * (X, X ) is the cofactor of the element. This identity allows us to remove four of the elements of S if we know 
the determinant \S\ and there remain two unitary conditions upon the magnitudes of the remaining, independent, 
elements. Thus in the end we see that 5, like U, is parameterized by nine real numbers; three magnitudes and six 
phases though the phase of the determinant may be stationary if the Hamiltonian in that basis is traceless. 



A. The two flavor approximation 

We have, so far, described everything in terms of three flavors but a quick scan through the literature by the reader 
will reveal that many studies have used only two. In this section we show how the three flavor S matrix formalism 
using the adiabatic basis can be separated in pieces which contain only two flavors. 

The motivation for the reduction in number of flavors one must consider comes from observation of the structure 
of the adiabatic Hamiltonian, equation (48). If we make the assumption that only one is significant at any given 
location, never two (or all three) simultaneously, then neutrino mixing occurs only between states \vi) and \vj) and 
the third state, \i>k), is decoupled. It is not immediately obvious that only one is significant at any given location: 
each Tij is proportional to the same derivatives of the potential dd V a fj / dx and the only difference is the matter mixing 
angle prefactors and the difference between the eigenvalues Skij in the denominators. Nevertheless, this is often the 
case for most density profiles. 

We can express the full Hamiltonian H ^ as the sum of three terms 

H {a) = #(i2) + #(13) + #(23) (62a) 

r 12 e 2l7TSct > 12 0\ si /0 Tige 2 "^ 13 



r^e" 2 "^ 13 




(62b) 

and then introduce the three S'-matrices £(12)) $(13) an d £(23) we obtain from the substitution of each into 
equation (56). Due to the structure of each the three S^j) have the form 

/ C(12) 77(12) \ 

S(i2) = -?7(i2) C* 12) (63a) 
\0 \j 

C(13) 7?(i3) 

5(i3) = ( 10^ j (63b) 

-J 7(13) C(*13) 

1 

S{2Z) = ( C(23) 77(23) I (63c) 



77(23) C(* 



(23) 



where the Qij) 's and 's are Cayley Klein parameters. Now it is simply a case of adapting this approximation to 
the situation at hand. In figure (1) we see that the L resonance always involves mixing between states \i>\) and \v2) 
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so if the density profile under consideration possesses only an L resonance then the S matrix describing the evolution 
of the neutrinos through the profile will have the structure of £(12) • The H resonance mixes states |i> 2 ) and l^) for 
a normal hierarchy so for a profile containing just an H resonance the S matrix will have the structure of £(23) • By 
applying the same reasoning for all the different possibilities we can assign the appropriate for any resonance 
shown in the figure. For a density profile possessing multiple resonances we can apply the group product rule for 
the evolution operator, equation (59), to divide the profile into sub-domains such that, within each, there is just one 
resonance. Since the structure of the 5* matrix for each sub-domain is given by the reasoning above the S matrix for 
the entire profile is then the time-ordered product of the appropriate <%j)'s. 

Once we have decided which pair of states are mixing we can then either solve the reduced problem 

t dsW =s m s{ 2) (64) 

dx to) V ' 

where S^ 2 ' is a 2x2 matrix and 

using the expression for given in equations (50a)- (50c), use a straight two flavor calculation a la Kneller and 
McLaughlin [16], or utilize some other alternative or approximate method. However one determines S^ 2 \ once it has 
been found one simply constructs the appropriate three flavor S matrix (or matrices in the case profiles with L, H 
and/or fir resonances) as in, for example, Kneller, McLaughlin & Brockman [19]. The disadvantage of applying a 
series of two flavor approximations is that some phase information can be lost. Nevertheless, many features of a flavor 
transformed neutrino signal can often be determined in this way. 



B. Identities of the Scattering Matrix 

The problem of neutrino propagation through supernovae has received considerable attention the in the past few 
years. The signal from the next Galactic supernova has the potential to reveal a great deal of information about 
both the supernova and the mixing parameters for the neutrinos. For example, if the angle #13 is not too small then 
dynamic MSW effects may be observed and some authors have also considered the possibility of observing effects from 
a non-zero CP phase, ft is upon this possibility of observing the non-zero CP phase that we now focus our attention. 

In the adiabatic and matter basis we have shown that the Hamiltonian is independent of the CP phase e if = 
V T = 0. This occurs because the eigenvalues and 6*i 2 and #13 are all independent of e in this limit so that the 
non-adiabaticity parameters T^, Ti3 and r 23 are also independent of e. If the Hamiltonian is independent of e 
then the S matrix must also be independent of e and, therefore, we must have S^(X, X , e) = S^(X,X ,0) and 
SH(I,I 0l e) = S^(X,X ,0). As a result, all the survival and crossing probabilities for the matter or adiabatic 
neutrino states, P{\vj) — > \vi)) = |<S^" i,a ^| 2 , are also independent of e. Any dependence upon e for the survival/crossing 
probabilities of states in other bases can only enter explicitly in the unitary transformation to those bases. 

For example, let us consider the transformation to the flavor basis. The flavor basis survival/crossing probabilities 
are found from the matrix related to 5< m )(A, X ,e) and S< a )(A, X ,e) by 

SW (X, X , e) = U(X, e) (X, X , e) W(X , e) = U(X, e) W(X, e) S^(X, X , e)W\X , e)W(X Q , e). (66) 

In above expression we can set W^(Xq, e) = 1 if all the phases in equation (46) are set to zero at Xq. However when 
splitting the S matrix into pieces, as in equation (59), then one should retain this phase information, so = 1 only 
for the first, rightmost, S matrix. 

Since S (m) (X, X Q , e) = 5 (m) (A, A o ,0) the following identity must be obeyed: 

U\X, e) S U) (X, X , e) U(X , e) = W{X, 0) S (J) (X, X , 0) U(X , e = 0). (67) 

Now let us consider how the CP phase enters into U. We saw from equation (18b) that U could be written as 
U = B0, 5) ©23(^23) E{e) 613(^13) ©12(^12) A(6ii, 012, 0.3) but, when = V T = 0, we have also seen that both 9i 2 
and #13 are independent of e and, furthermore, we have also made the decision that di, a 2 and 0:3 are constants. This 
means ©12(2;, e) = @i 2 (a;,0), ©13(2;, e) = ©13(2;, 0) and A(x, e) = A(x,0) so when we insert the decomposition of U 
into equation (67) we find that what survives can be written as 

B\X,e) SW(X,X ,e)B(X ,e) 

= © 23 (A, e) E(X, e) © 23 (A, Q)B\X, 0) S«\X, X , 0) B(X , 0)© 23 (Ao, 0) tf(X , e) ©^ 3 (A , e). (68) 
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The two combinations 9 23 (A, e) E(X, e) e\ 3 (X, 0) and e 23 (A , 0) E\X ,e) ®l 3 (X , e) arc 

e 23 (A,e)£(A,e)e 23 (A,0) = 

1 

S 23 (A, e)§23{X, 0) + c 23 (A, e)c 23 (A, 0)e"W S 23 (A, e)2 23 (A, 0) - c 23 (A, e)S 23 (A, 0)e^ x ) | , (69) 

S 23 (A, e)S 23 (X, 0) - S 23 (A, e)c 23 (X, 0)e"W c 23 (A, e)£ 23 (A, 0) + S 23 (A, e)S 23 (A, 0)e l? « 

e 23 (X ,0)St(X ,e) et 3 (x ,e) = 

1 

S 23 (A ,e)S 23 (A ,0) +£ 23 (X ,e)5 23 (Xo,0)e^ x «) 5 23 (A , e)S 23 (A , 0) - S 23 (A , e)2 23 (X , 0)e^ ? ™ 

S 23 (X ,e)2 23 (X ,0) - 5 23 (X ,e)S 23 (X ,0)e-"( x «) c 23 (A , e)2 23 (A , 0) + S 23 (A , e)S 23 (A , 0)e~^ x «) 

(70) 

respectively. After inserting these matrices into equation (68) we find nine relationships between the elements of 

S(f\X,Xo,e) and the elements of S^(X,Xo,0). The flavour basis survival/crossing probabilities P a fj = \S^\ 2 

may be found for both S^\X, A , e) and (X, A , 0) and using the relationships between the elements of the two 
matrices we find that we can derive four, non-trivial identities for the probabilities: 

• Pee(X,X ,e) = P ee (X, A ,0), 

. P ef ,(X, X ,e) + P eT (X, X ,e) = P eM (A, X , 0) + P er (A, A , 0), 

• P^ e (X, X ,e) + P Te (X, X ,e) = P^X, A , 0) + pH ] (A, A , 0), 

• P„„ (X, X , e) + P^ (X, X , e) + P Tfl (X, X , e) + P TT (X, X , e) = P„„ (X, X , 0) + P^ (X, X , 0) + P T/1 (X, X , 0) + 
P TT (A,A o ,0). 

The first two identities were also found by Balantekin, Gava and Volpe [56]. 

However when computing the fluxes at Earth one requires a slightly different set of probabilities: the probability 
that an initial flavor state emerges as a given matter/mass state. These probabilities are found from the matrix S^ m ^ 
given by 

5 (m/) {X, X ,e) = S (m) {X, X Q , e) (}\X Q , e) = W(X, e) S {a) {X, X , e) W*(X , e)U*(X , e). 

(71) 

Again the invariancc of (A, A , e) with regard to e when V„ = V T = means that we have the following identity: 

S^(X, A , e) L>(A , e) = S^(X, A , 0) U(X , 0). (72) 

which ultimately leads to 

(A, A , e) B(X , e) = S^ mf \X, X , 0) B(X , 0) 9 23 (A , 0) E^X , e) 9 23 (A , e). (73) 

This is very similar to equation (68). Again we have nine relationships between (A, A , e) and S^ mf \X,X o ,0) 

which lead to the following, non-trivial identities for the probabilities P ia = \S^^\ 2 : 

. P le (A,A ,e) = P le (A,Ao,0), 

. P 2e (A,A ,e) = P 2e (A,A ,0), 

. P 3e (A,A ,e) = P 3e (A,A ,0), 

. P 1(U ( A, A , e) + P lr ( A, A , e) = P 1(U ( A, A , 0) + P lr ( A, A , 0) , 

• P 2li ( A, A , e) + P 2T (X, X , e) = P 2[i (X, A , 0) + P 2r (X, X , 0) , 
. P 3AJ ( A, A , e) + P 3T ( A, X , e) = P 3ai ( A, A , 0) + P 3r (X, A , 0) . 
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At the detector on Earth the supernova neutrinos are measured as flavor eigenstates thus, for example, the electron 
neutrino signal is the appropriate linear combination of survival probabilities multiplied but the original fluxes. One 
can see from these identities above that the CP phase effects will not show up in observed supernova neutrino signal 
if the u^jVt fluxes emitted from the supernova neutrino sphere are equal [55, 56]. However, as pointed out in [56] if 
the Up, v T fluxes are not equal when they are emitted, then there is no guarantee that the fluxes as observed in the 
detector will be independent of the CP phase. 

In the case where when e ^ and 5V^ T =/= these identities no longer apply because the eigenvalues and the 
mixing angles now become functions of e. One might imagine that some sort of cancellation occurs such that the 
three non-adiabaticity parameters are independent of e. However an examination of the adiabatic basis shows that a 
complete cancellation cannot occur. For all three non-adiabaticity parameters, equations (50a)-(50c), there emerges an 
imaginary component proportional to sin e dSV^r / dx and this term cannot be canceled by concomitant changes in the 
eigenvalues and/or mixing angles. For all three Tij we have a situation where T^e) ^ 1^(0) thus H^(e) ^ H^(0), 
S^(e) ^ 5^(0) and, finally, P^e) ^ P i3 (0) for any basis. 

VI. SUMMARY AND CONCLUSIONS 

We have considered the problem of a generalized, 3-flavor, neutrino mixing that includes matters potentials for both 
ll and t flavors and CP violation. We presented expressions for the eigenvalues and matter mixing angles and pointed 
out that mixing phases that were zeroed in the vacuum are not necessarily zero in matter. We found that in the limit 
that the mu and tau potentials are equal, the eigenvalues and matter mixing angles #i 2 and #13 are independent of 
the CP phase, although 623 does have CP dependence. 

We introduced the 3-flavor adiabatic basis. In this basis the Hamiltonian is completely off-diagonal and the behavior 
of the neutrinos is largely determined by the 3-flavor adiabaticity parameters. This is a useful basis in which to 
calculate neutrino flavor transformation; further it gives a straightforward picture of the effects of the CP phase. In 
the limit that the mu and tau neutrino potentials are the same, the Hamiltonian in the adiabatic basis is independent 
of both 623 and the CP violating phase, so CP phase effects appear only in rotations into and out of the flavor basis. 

We discussed the 5* matrix for three flavors, as well as two flavor S-matrix approximation. Using the three flavor 
S-matrix, we found several non-trivial identities related to the observability of the CP phase. In a future study [66] 
we shall discuss how one can formulate efficient algorithms for the calculation of S and present some calculations of 
three flavor oscillations made with them. 
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